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Abstract 

We prove that the twist-3 chiral-odd parton distributions obey simple Gribov-Lipatov- 
Altarelli-Parisi evolution equations in the limit N c — > oo and give analytic results for the 
corresponding anomalous dimensions. To this end we introduce an evolution equation for 
the corresponding three-particle twist-3 parton correlation functions and find an exact an- 
alytic solution. For large values of n (operator dimension) we are further able to collect all 
corrections subleading in N c , so our final results are valid to 0(1/N 2 ■ ln(n)/n) accuracy. 



PACS numbers: 11.15.Pg, ll.10.Hi, 12.38.Bx, 13.85.Qk 



1 



The increasing precision of experimental data from LEP, HERA and the TEVATRON 
requires understanding of higher twist corrections induced by correlations of partons in 
the colliding (produced) hadrons. The twist-3 parton distributions play a distinguished 
role in spin physics, where they can be measured as leading effects responsible for certain 
asymmetries. In particular, the E143 collaboration has reported the measurement of the 
twist-3 contribution of the polarized structure function g2(x,Q 2 ) [|J, thus providing with 
the first experimental test of quark-gluon correlations in the nucleon. The chiral-odd parton 
distribution tiL{x,Q 2 ) is expected to be measurable in the Drell-Yan production @||. 

The Q 2 evolution of twist-3 distributions is usually believed to be quite sophisticated due 
to mixing with quark-antiquark-gluon operators, the number of which increases with spin 
(moment of the structure function). For the case of the chiral-even twist-3 flavor-nonsinglet 
distribution g 2 (x,Q 2 ) it has been observed by Ali-Braun-Hiller |4j that the operators in- 
volving gluons effectively decouple from the evolution equation in the large N c limit, that is 
neglecting 0(1/N 2 ) corrections ||. Consequently, g 2 (x,Q 2 ) obeys a simple GLAP evolution 
equation and the corresponding anomalous dimension is known in analytic form ||. The 
statement holds true with full account for effects subleading in N c but for large moments n. 
Thus the claimed accuracy is in fact 0(1/N 2 ■ ln(n)/n). The ABH evolution equation gives 
a guide to the expected small x and large x behavior which is important for experimental 
extrapolations, and is used to rescale the model predictions to high values of Q 2 of the 
actual experiments. 

Physics of this decoupling is so far not understood, and it is probably not related to 
usual simplifications of the iV c — > oo limit. The observed phenomenon appears, however, to 
be quite general. In this letter we demonstrate that the same pattern is obeyed by chiral- 
odd parton distributions h L (x,Q 2 ) and e(x,Q 2 ) as well, albeit with different anomalous 
dimensions. For all practical purposes this solves the problem of the Q 2 evolution of twist-3 
nonsinglet parton distributions, since the corrections 1/N 2 are small. 

Following |3],[7| we define the parton distributions in question as nucleon matrix elements 
of nonlocal light-cone operators §| 
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(PS\i){z)i){-z)\PS) = 2M J^dxe 2l(Pz)x e(x,fi 2 ) (1) 

and 

(PS\4>(z)a^z u i l5 ij(-z)\PS) = ^S^(P-z) J^dxe 2i ^ x h x {x^ 2 ) 



IMz,^- J^dxe 2i ^ x h L (x,^). (2) 



(P 

Here z is a light-like vector z 2 = 0, P and 5* are the nucleon momentum and spin vectors 
(p2 = M 2^ S 2 = _ M 2 ; P .s = Q) t respectively, and S ± ^ = S„- P^S ■ z)/(P ■ z) + M 2 Zfl (S ■ 
z)/{P ■ z) 2 . A thorough discussion of the parton interpretation and of physical relevance of 
these distributions can be found in Ref. |§. 

As it is well known, equations of motion allow to express the twist-3 quark- ant iquark 
distributions in terms of quark-gluon correlations. This is visualized most explicitly in form 
of the operator identities (see Eqs.(27),(28) in ||) 

$(z)ij)(-z) = tp{0)ip{0) + f 1 du f U dt S(u, t, -u) , (3) 

JO J~u 

^(z)a^z„i^(-z) = [vp(z)a^z u i^(-z)] tw i st 2 + iz^ udu / tdt S(u,t, -u), (4) 

JO J—u 

where we have introduced shorthand notations for the nonlocal operators 

S(u,t,v;/j, 2 ) = i)(uz)a IJi gG vi (tz)z v ,z v ip{vz) , 

S(u, t, v; fj 2 ) = tjj(uz)ia IM ^ 5 gG^(tz)z IM z u i)(vz) . (5) 

The relations in (||) and (f|) are exact up to twist-4 corrections and neglecting operators 
containing total derivatives §J which are irrelevant for our present purposes. 

Note contribution of the local scalar operator in Eq. (|3]). It gives rise to the nucleon o- 
term and does not have a partonic interpretation. The remaining contributions give rise to 
three-particle parton distributions describing interference between scattering from a coherent 
quark-gluon pair and from a single quark. In particular, one can define 

(PS\S(u,t,-u)\PS) = -8M(p-z) 2 [da [ d(3e i{p - z)[a{u ~ t)+ ^ u+t)] D g (a,(3), (6) 



and a similar quantity D g (a, (3) for S. The support properties of D g (a, /3),D g (a, j3) and their 
parton interpretation are discussed in Ref. [[II]]. The variables a, (3 and a — (3 have physical 
meaning of the momentum fractions carried by the antiquark-, quark- and gluon-partons, 
respectively, and D g ,D g vanish unless \a\ < 1, |/3| < 1, and \a — (3\ < 1. Combining (||), 
and @ one can express the twist-3 structure functions h^x) and e(x) in terms of a certain 
integral of these quark-antiquark-gluon distributions |13| . Note that only such integral (over 
the gluon momentum) is potentially measurable in inclusive reactions like DIS or Drell-Yan 
processes. 

The Q 2 dependence of the twist-3 distributions is governed by the renormalization group 
(RG) equation for the corresponding nonlocal operators S, S. To leading logarithmic accu- 
racy the evolution of S and S is the same; hence we drop the "tilde" in what follows. 

We find it convenient to use a general approach of |T| to write RG equations directly 
for the nonlocal operators. To this end we introduce the Mellin transformed operators fllH 



S(u, t, v) 



1/2+ioo 



dj(u-vy- 2 s( 3 ,o, £ 



u + v -2t 



u — v 



(7) 



2m Ji/2-i 

Here j is the complex angular momentum; operators with different j do not mix with each 



other. The Mellin transformed operators satisfy the RG equation [|Tl 



' d d\ 



a s 
2^ 



(8) 



where the kernel Kj(^,rj) is subject to an explicit calculation. Neglecting all contributions 
which are down by 1/N 2 we get 



1 \ 5 a^ ^ 6>(e - 77) 1 + 77 



N, 



5 drj! 1 + 77/ 



-ie - 1 + i 



1 + 



1 — T] 



+ 



1 - rj/ / 



1 drjt l — 



1 + 



1 + 77/ 



-4 



-4 



(i + O 3 

oijj - 
(i-O 3 



i-O 
1 + 77J 



1—7] 

1 + 77 
J + l 



.1 + 77, 



(9) 



The kernel turns out to be very similar to the corresponding kernel for the evolution of 



chiral-even twist-3 operators in Ref. M. Derivation of (pi) will be given elsewhere |13 



To solve the RG equation in (Rf) we consider the conjugate homogeneous equation 



(10) 



Each eigenfunction of (S) corresponds to a multiplicatively renormalizable nonlocal operator 
with the anomalous dimension 7^: 

rl 



(11) 



where 6 = (lliV c — 2Nf)/3. To prove this, multiply (§) by (/>j(£) and integrate over £. Using 
(|T0|) one obtains the conventional diagonalized RG equation with the solution in flTT|). 
We were able to find two solutions for the equation in flT0| ) with eigenfunctions 



+ (O = i, <T(0 = £ 



(12) 



(they do not depend on j, so we drop the subscript). The corresponding eigenvalues (anoma- 
lous dimensions) equal 



7 } = 2N.^ + l)+y E -\-^ j 

77 = av«{*0 + i) + ^-j + 5oTT) 



(13) 
(14) 



where ?/>(z) = 4- lnT(2;) and 7^ is the Euler constant. Validity of Eqs. ( |i3|) and ( |I~4"D can be 



checked by a straightforward calculation. By comparison with results of the numerical eval- 
uation of anomalous dimensions for integer j = n (see below) we conclude that our solutions 
always correspond to operators with the lowest anomalous dimension in the spectrum. 

The superscript ± corresponds to the "parity" under £ — > — £: due to the symmetry of 
the kernel one can look for separate solutions in the space of functions which are even (odd) 



with respect to the substitution £ — > — £. From Eqs.(|3|) and (|), one sees that the relevant 
quantities for e(x) and hi(x) are even and odd "^-parity" pieces of the nonlocal operators. 

Substituting the definition (|7|) into (|j) and (f|), and changing the integration variable 
t to £ = — t/u, we observe that the nonlocal operators in (|IID with the particular choices 
of the weight function (|T2p are precisely those which give rise to twist-3 quark- ant iquark 
operators at the tree level. Taking the nucleon matrix elements, we get for the moments, 

M n [e](Q) = L^' b M n [e](fi), (15) 
M n [h L ]{Q) = L^/ b M n [h L ](fx), (16) 

where M. n \hL\ = j\dxx n JiL{x), A4 n [e] = j\dxx n e(x), and L = a s (Q)/a s (fi). Hl is the 
genuine twist-3 contribution to h,L, after subtracting out the twist-2 piece ||. The solutions 
for the evolution in ([15]), ( ]IBD present the principal result of our paper. 

Expansion of nonlocal operators at small quark-antiquark separations generates the series 
of local operators of increasing dimension. In particular, expansion of S(j,£) for positive 
integers j = n > 2 generates the local operators 9 n k = ^p(D ■z) k ^ 2 aG{z- D) n ~ k ?p (k = 
2, 3, ■ • ■). The substitution £ — > — £ for the nonlocal operators corresponds to k —>■ n — k + 2 
for the local operators. From the kernel @ one can calculate the mixing matrix for the 
local operators even as well as odd under k — > n — k + 2. One can check that the result for 



the odd case coincides with Eqs. (3.14)-(3.16) in Ref. [12| to the stated accuracy. In 

particular, neglecting the 1/N% terms in the mixing matrix of [0 we obtain the evolution 
equation (|16| ) with 



- ! / 2k 



Mn[h L ](n) = E 1 - —To) KM- (17) 
where 6 n ,fc(/^) {k — 2, ■ ■ ■ , [(n + l)/2]) are reduced matrix elements of the independent quark- 



ant iquark-gluon operators in the notation of [§,[T3I • This is the consequence of the fact that 
in this limit we have 



1 1 ' 2k \ .. / 21 



Sl'-^r"" 1 -^'* (18) 



k=2 



where Xki is the mixing matrix in the notation of [12 as 



r 21+11 

f d d\ a [ 2 J 

+ m d~g) KM = 27 £ ^M/^)- (19) 

The coefficients (1 — 2k/ (n + 2)) correspond to of (0) in the nonlocal operator lan- 

guage. Equation flTTD is precisely the operator giving Hl at the tree level. All operators with 
higher anomalous dimensions decouple from the Q 2 evolution in the iV c — > oo limit. 

The complete spectrum of anomalous dimensions in the N c —>■ oc limit obtained by the 
numerical diagonalization of the mixing matrix in |12| is shown in Fig. 1, together with our 
analytic solution for the lowest eigenvalue. 

To illustrate numerical accuracy of the leading-iV c approximation, consider the exact 
result [|l2j (including 1/AT? corrections) for the evolution of the n = 5 moment of Hl, which 
is the lowest moment in which mixing appears: 

M 5 [h L ](Q) = [0.4166 5 , 2 (/i) + 0.1936 5 ,3^)]^ 12 - 91/fe 

+ [OmSbsM - 0.0506 5 , 3 (/x)] L lsm ' b . (20) 



This is reduced in the large iV c limit to 



M 5 [h L ](Q) 



3 1 

^&5,2(m) + 



L 13.7/ fe _ (21) 



One observes: (i) the lowest anomalous dimension and the coefficients in front of the two 
terms in this contribution are close to their exact values; (ii) the contribution of the operator 
with the higher anomalous dimension is small (~ 1/N 2 ). It is the latter observation which is 
crucial for the phenomenological importance of our results, since it means that description 
of each moment of the twist-3 distribution requires one single nonperturbative parameter. 
A similar check for e(x) [i.e., (18)-(21)] has also been done in [[14]]. 

We can make this point even stronger, by observing that admixture of operators with 
higher anomalous dimensions is suppressed at large n for arbitrary values of iV c . The full 
evolution equation with account of all 1/N 2 terms is complicated and will be given elsewhere 
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T3fl . However, it simplifies drastically in the limit j ^ oo and coincides with the large- j 



evolution equation considered in [|J. We get 

iM) = hi + iE - 3/4] m + AW({) 

+ Ne (izlV m - m + Nc i\ fi+iV «t) - m , (22) 



-l \i — t] J £, — v + v J v ~ £ 

with Cf = (N 2 — l)/2N c . Comparing to the large- j limit of the kernel in (|9]), the only 
difference is in the replacement 2N C — > 4C/ in the first term. Thus, the functions in flT2|) 
still provide the solution, and the anomalous dimensions are shifted by 

lf^lf + (4C7 - 2AQ [In j + 1E - 3/4]. (23) 

With this modification of the anomalous dimensions, the results in flip]) are valid to the 



0(1/N^ ■ ln(n)/n) accuracy. 

To summarize, solutions in the present paper provide a powerful framework both in 
confronting with experimental data and for the model-building. From a general point of 
view, they are interesting as providing with an example of an interacting three-particle 
system in which one can find an exact energy of the lowest state. For phenomenology, 
main lesson is that inclusive measurements of twist-three distributions are complete (to our 
accuracy) in the sense that knowledge of the distribution at one value of Q% is enough to 
predict its value at arbitrary Q 2 , in the spirit of GLAP evolution equation. This allows to 
relate results of different experiments to each other, and to compare with model calculations 
which typically are given at a very low scale. 

As shown above, the 1/N 2 corrections are not large for n = 5 and further decrease as 
ln(n)/n at large n. Thus, the only possibility for sizeable 1/N% effects might be for small x 
behavior, in case the location of the singularity in the complex j plane of the exact anomalous 
dimension happens to be on the right of the singularity of the leading iV c result. A detailed 
study of the small x behavior goes beyond the tasks of this letter. 

We thank the Institute for Nuclear Theory at the University of Washington for its hos- 
pitality and the DOE for partial support during our visit which initiated this study. 
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FIGURES 




FIG. 1. Spectrum of anomalous dimensions of twist 3 operators for h^x) in the limit N c — ► oo. 
Fhe solid line shows the analytic solution (14). 
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